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AUTOMATED FAULT DIAGNOSTICS FOR FREQUENCY- 
DEPENDENT CIRCUITS BY CORRELATION AND VARIANCE TECHNIQUES 

By Henry B. Metcalf 
Electronics Research Center 
Cambridge , Massachusetts 


SUMMARY 

A computer program, NASAP (Network Analysis for Systems 
Applications Program) , has been adapted to perform fault diagnosis 
in frequency dependent systems . 

The program simulates a short circuit and an open circuit 
for each network parameter and prints out a modified transfer 
function of the network corresponding to each simulated failure. 
Experimental and simulated data are subjected to correlation and 
variance analysis in order to rank all parameters according to 
probability of failure in the network model. 

The program lists, at present, single failures in order to 
decreasing probability. Possible extensions of the program to 
multiple failures are examined. 

I.- INTRODUCTION 
Review of Work 

The computer program, NASAP, has been adapted to perform 
fault diagnosis in frequency-dependent circuits. The present 
version of NASAP can handle linear circuits, which consist of 
constant-value passive elements, and independent or dependent 
current and voltage sources. The program has been written in 
CDC FORTRAN IV to make it easy for the intended user to under- 
stand the algorithms and to make changes within the program. The 
version of NASAP described within this report has been run on the 
IBM 7094. 

The technique employed in the failure diagnostic method has 
evolved from three sources. The first of these is linear-graph 
theory and flow-graph theory. Each circuit is given a linear 
graph description which is then formulated into flow-graph terms 
by a computer subroutine thus forming the basis of the computer 
algorithms for later output requests. 

The second area is that of correllation and variance analysis. 
The correllation coefficient is the statistical index used to 
measure the similarity observed between like sets of data. It 
is used here to provide a standardized measure of the degree of 



deviations observed in the simulated failures. The variance 
measure is employed to provide an index of the percentage devia- 
tion caused by a simulated failure. 

The last source and also the foundation of the technique is 
NASAP . The fault-isolation technique developed in this report 
has been formulated for use with NASAP. 

Strategy for Identification of Faulty Elements 

For a given network, a computer program simulates the failure 
of a component in a multi-element network. This simulation is 
accomplished by modifying NASAP as follows: 


1) Solve for the transfer function of the network of interest, 
designated by F(0). 

2) Compute F (k) for 2N distinct cases; namely, when each 
of the N elements is short-circuited and open-circuited. 

The functions F (k) with k = 1 , 2 , . . . , 2N but k ^ 0 define 
a set of modified transfer functions where F(2J-1) are 
short circuits and F(2J) are open circuits, where 

J = 1 ,2 , . . . ,N. 

3) Evaluate F(0) and F(k) numerically at specified complex 
frequencies by taking the absolute value of the transfer 
functions and the modified transfer functions; this 
yields 2N sets of data for the modified transfer func- 
tions, F (k) , and one set of data for the original trans- 
fer function, F(0). 


4) Solve for the correlation coefficients and the variance 
coefficients between the original transfer function and 
the modified transfer functions. This process yields a 
table of identification coefficients for all possible 
failures . 

5) Measure the correlation and variance coefficients between 
the observed response and the desired response, when a 
variant response is observed in a circuit. These coeffi- 
cients are then compared to the table of identification 
coefficients and the possible failures are printed out 

in their order of most probable occurrence. 

Computational Procedures for Locating Fault Elements 
To pinpoint the failure from a knowledge of: 
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F ( 0 ) Transfer Function 
F (k) Modified Transfer Functions 



two coefficients are derived. 

a) Correlate F(0) with each of the F(k) modified transfer 
functions to yield the correlation coefficient. 

b) Measure the variance caused by each short circuit or 

open circuit at the specified frequencies and the variance 
coefficient results. 

c) Both coefficients can be calculated for every F (k) . 

Experimentally the identification is performed as follows: 

1) Construct the network experimentally and measure the 
transfer function. 

2) Determine how the modified transfer functions deviate 
from the calculated transfer function. 

3) Use criteria (a) or (b) or both to identify the faulty 
element . 

4) List a small number of elements that are expected failures 
in their order of most probable occurrence. 

Application of Results 

The computer program output provides the following: 

1) The transfer function of the network of interest. 

2) The modified transfer functions corresponding to solu- 
tions for each short circuit and open circuit. 

3) The numerical evaluation of each transfer function at 
specified frequencies. 

4) The correlation coefficient between the calculated trans- 
fer function and each modified transfer function. 

5) The variance coefficient computed for each modified 
transfer function. 

6) The most probable failures in terms of the experimental 
results . 

7) The program is written in CDC FORTRAN IV and is presently 
being run on the IBM 7094. The program is available from 
project COSMIC. 
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II.- METHODS 
The Topology Equation 


The basic algorithm employed in the computer program is the 
topology equation which provides a method for finding the desired 
cause-effect relationship between an input and output variable. 

All applications of the topology equation in this report deal with 
closed-signal flowgraphs. 

The closing element P is the performance parameter. Q 
refers to the element where failure is being simulated and the 
topology equation, H = 0, can be expanded as follows: 


H(P,Q) = H(P,Q) + PH (P ' , Q) + QH (P , Q 1 ) + PQH(P',Q') 


It states that P can be expressed in terms of Q, since 
H(P,Q) = 0 is linear in P and Q. Although the topology equation 

applies to all types of linear systems, it can be illustrated by 

a flowgraph (such as the one) in Figure 1. 

In the case of flowgraphs the topology equation yields: 

H (P ' ) - loops with the P tag present 

H(P) - loops that are deprived of the_P tag 

and similarly for H(Q') and H (Q) . 

The combination of loops in the topology equation of a 
closed network is given by: 

H = (-l) m H(m) where 

H(m) = L(m,n) and 

m is the order of the loop and n is the number of loops belonging 
to each order. L(m,n) is the nth loop of order m. A loop is 
defined as a sequence of transmittances such that every node is 
common to two and only two transmittances of the loop, one termina- 
ting at the node and one emanating from the node. In a first- 
order loop every node can be reached from every other node. A 
loop of order m is a set of N disjoint first-order loops. 

For Figure 1 the loops are listed in Table I. The topology 
equation then is: 

H(P,Q) = 1 - MN 
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H (P ' , Q) - -PKL 
H (P , Q ' ) = JQR 
H (P' ,Q' ) = JKL QPR 



Figure 1.- Loop example 
TABLE I.- LISTING OF LOOPS FROM FIGURE 1 



The Modified Transfer Function 
Solving H = Q for 1/P gives 

ik = kL kE—2 . = f ( o ) , the transfer function 
H (P ) 

H (P ' ) = H(P' ,Q) + QH(P' ,Q ' ) 

H(P) = H(P,Q) + QH (P , Q ' ) 
implies 

1 _ H (P ' ,Q) + QH (P 1 , Q 1 ) = p(0) 

P H(P,Q) + QH (P , Q 1 ) 
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Table II defines the variance in impedance short circuit 
and open circuit of parameter Q, as derived from the above 
formula for F(0). 


TABLE II.- MODIFIED TRANSFER FUNCTIONS 


Cause of 
Variance 

Q = Y 

Q = z 

Short 

Circuit (SC) 

1/Q = 0 

o 

II 

o 

Open 

Circuit (OC) 

o 

II 

a 

1/Q = 0 

Y° C = Z SC 

H (P 1 ,Q) 
H(P,Q) 

H (P 1 , Q) 
H(P,Q) 

Y sc = z oc 

H (P 1 , Q ' ) 
H (P , Q ' ) 

H (P 1 , Q 1 ) 
H (P , Q ' ) 


Illustrative Example for Numerical Data 

Experimental data of the transfer function F(0,i) and the 
modified transfer function F(k,i) at specified frequencies w(i) 
are printed out by the computer program. This will be illustrated 
for the circuit in Figure 2. 



Figure 2.- Illustrative circuit 
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The circuit was constructed with components listed in 
Table III and I g/Vyi was measured experimentally. 

TABLE III.- COMPONENT VALUES FOR ILLUSTRATIVE NETWORK 


ELEMENT 

VALUE UNIT 

VI 

1.000 

volt 

Rl 

1.000 

K ohm 

L7 

5.412 

MH 

R7 

1.262 

k ohm 

C7 

200 . 000 

PF 

L6 

5.645 

MH 

R6 

1.307 

K ohm 

C6 

212.800 

PF 

R3 

1.228 

K ohm 

C2 

281.800 

PF 

R5 

1.254 

K ohm 

R8 

119 . 000 

ohm 


Table IV gives the data sets for F (0) , the original trans- 
fer function, and F(l), F(3), F(4), three modified transfer 
functions. F ( 1 ) corresponds to the short circuit of Rl , F(2) 
corresponds to the open circuit of Rl , and F(3) is the short 
circuit of L7 . 

Since F(2,i) = 0, the data set of Rl open circuit, it is 
not listed in the table. 

It is important to note that in many of the following tables 
the data have been carried to five or six places. This is for 
illustrative purposes only since the input data imply only four- 
figure accuracy at the most. 


TABLE IV.- DATA SETS 


w (i) 

F (0, 

i) 

F ( 1 , 

i) 

F (3, 

i) 

F (4, 

i) 

10 1 * 5 

. 362874 

E-03 

.193424 

E-03 

.156740 

E-03 

.632455 

E-08 

io 2 

.362874 

E-03 

.193424 

E-03 

.156740 

E-03 

.200000 

E-07 

10^.5 

. 362873 

E-03 

.193424 

E-03 

.156740 

E-03 

. 632455 

E-07 

io 3 

. 362879 

E-03 

.193423 

E-03 

.156740 

E-03 

.200000 

E-06 

10 3 ' 5 

.362824 

E-03 

. 193419 

E-03 

.156739 

E-03 

.632455 

E-06 

io 4 

.362378 

E-03 

.193376 

E-03 

.156630 

E-03 

.199998 

E- 05 

10 4 .5 

.357991 

E-03 

. 192952 

E-03 

. 156641 

E-03 

.632382 

E-05 

io 5 

. 320148 

E-03 

.188740 

E-03 

.155739 

E-03 

.199782 

E-04 

10 s * 5 

.160346 

E-03 

.150203 

E-03 

.145178 

E-03 

.629477 

E-04 

io 6 

. 124238 

E-04 

. 204408 

E-04 

.387241 

E-04 

. 342242 

E-04 
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The Correlation Coefficient 


The correlation coefficient (R) for F(0,i) and some set 
F(k,i) is the mean of the paired products of the deviations of 
each score from their respective means when these deviations are 
measured in units of their respective standard deviations, 

D (0 ,i) and D (k,i) . 

F(0,i) - G ( 0 , i ) Vf ( k , i ) - G(k,i) \ 

, D(0,i) A D(k,i) / 

F (0 ,i) F (k , i ) - G ( 0 , i ) G(k,i) 

D(0,i) D (k , i) 



where the summation extends over the frequency interval specified 
by M measurements; hence, the means are: 


G(0,i) = Z F(0,i)/M, and similarly G(k,i) = Z F(k,i)/M 


and the standard deviation is : 

0(0,1) ^ F ^°’ i)2 - G(0,i) 2 )’ i 
Similarly , 

D (k , i ) F ^ k,1)2 " G(k ,i) 2 ) 35 for K ^ 0 


Each point evaluated corresponds to a point on a curve and 
each data set describes a curve; hence, a value of R is assigned 
for each short circuit and open circuit. That is, the deviations 
of each modified transfer function, F(k), from the original 
transfer function, F(0), result in a separate value of R as 
shown in Table V. The absolute value of R will ensure that it 
is positive. 
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TABLE V.- CORRELATION COEFFICIENTS 


k 

FAILURE 

SC 

oc 

R 

1 

Rl 

X 


.952979 

2 

R1 


X 

. 000 

3 

L7 

X 


.897051 

4 

L7 


X 

.794069 

5 

R7 

X 


.940948 

6 

R7 


X 

.794069 

7 

C7 

X 


.905976 

8 

C7 


X 

.938431 

9 

L6 

X 


.889213 

10 

L6 


X 

.741361 

11 

R6 

X 


.943269 

12 

R6 


X 

.741361 

13 

C6 

X 


.897955 

14 

C6 


X 

.929121 

15 

R3 

X 


.961076 

16 

R3 


X 

.000 

17 

C2 

X 


.000 

18 

C2 


X 

.948035 

19 

R5 

X 


.961537 

20 

R5 


X 

.000 


Variance Due to the Failure of a Single Component 

The correlation technique is not always sensitive enough 
to detect a failure. The variance coefficient, T, can be an 
alternative or is used in conjunction with the correlation 
coefficient, R. Both methods are measures of the deviation 
observed between the original transfer function and the modified 
transfer functions. When used together they make the fault 
isolation more effective. 

Table VI defines the variance, E, observed in P = 1/F(0) 
at the occurrence of a short circuit and open circuit of each 
element; that is, each F(K) . 

The variance defined in Table VI is computed for the example 
in Figure 2. The results for E(p,X) are listed in Table VII and 
the observed error is plotted in Figure 3 over the frequency 
range . 
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TABLE VI.- DEFINITION OF VARIANCE 


P (Q) 

H (P , Q) 

H (P 1 ,U) 

1 _ 1 1 
F (k) zsc yoc 

P (Q ' ) 

H (P , Q 1 ) 

1 _ 1 _ 1 

H (P ' ,Q ' ) 

F (k ) ZOC ysc 

E(p' , Q 1 ) 

mam 

1 1 
F (0) F (k) 

p 

1 

F (0) 

E (P 1 ,Q) 

P - P(Q) 
P 

1 1 
F ( 0 ) F (k) 

1 

F (0) 

E (P , Q ' ) 

1 1 
P P(Q') 

F (0) - F (k) 

1 

P 

F ( 0 ) 

E(P,Q) 

1 1 
P p¥ 

1 
p 

F (0) - F (k) 
F ( 0 ) 


From Table VI: 


F(0,i) - F (k , i) 
F(0,i) 


E{P,X) k 


where X is either Q or Q ' depending on whether we are dealing 
with a short circuit or an open circuit and k designates a 
particular modified transfer function. That is, if K = 1; 


E (P,x) k= i 


F(0,i) - F ( 1 , i ) 
F (0 ,i) 
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As before, F (0) is the original transfer function, F (1) is 
Rl short circuit, F(2) is Rl open circuit, F(3) is L7 short 
circuit, and F(4) is L7 open circuit. 

E (P , X ) k=2 = 1 over the frequency range, since F(2,i) = 0 
for all i and these data are left out of the table. 


TABLE VII.- FRACTIONAL VARIANCE OBSERVED AND 
VARIANCE COEFFICIENT 


u> (i) 

E (P , X) k = 1 



lO 1 ’ 5 

-.466968 

-.568060 

-.999983 

10 2 

-.466968 

-.568060 

-.999945 

io 2 ' 5 

-.466976 

-.568059 

-.999826 

io 3 

-.466976 

-.568066 

-.999449 

10 3,5 

-.466907 

-.568003 

-.998257 

o 
1 — i 

-.466368 

-.567496 

-.994481 

10 4 ‘ 5 

-.461015 

-.562444 

-.982335 

io 5 

-.410460 

-.513541 

-.937597 

5 5 

IqD.o 

-.063257 

-.094596 

-.607426 

io 6 

-.645297 

2.11693 

1.754730 

T(K) 

T(l) 

T (3 ) 

T(4) 


ZLog 10 |E(P,X) K=1 | 

ZLog 10 |E(P,X) K=3 | 

ILog 10 |E(P,X) K=4 | 


17.90326 

17.28895 

17.98819 
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Definition of the Variance Coefficient 


The variance coefficient, T, for the deviation between the 
original transfer function, F(0), and the modified transfer 
function F (k) is defined by B. 


D(P,X) k = 


F (0 , i) - F (k , i) 
F ( 0 , i ) 


(A) 


Where X is either Q or Q' , and k designates the modified 
transfer function being considered, and then: 


T = ELog 10 | E (P ,X) k | (B > 

where the summation is taken over all M in the frequency range. 
In Table VII for example, 

T(l> = ZLog 10 |E(P,X) k=1 | = ZLog 10 [ F( . 0 . ' . V ( -,i ) (1 . ^ i -I| 
and similarly for T(3) and T(4). Since: 


E (P ,X ) k _2 = If clearly T(2) = 0. 

Table VIII gives a complete listing of all variance coeffi- 
cients, (T) , and correlation coefficients, (R) , for the example 
in Figure 2. 


Fault Isolation Procedure 


In review: 

1) The network was constructed experimentally. 

2) The desired transfer function, F(0), was solved for. 

3) The modified transfer functions, F (k) , were established 
for all possible single failures. 

4) The correlation coefficients (R) between each modified 
transfer function and the desired original transfer 
function were formulated. 

5) Similarly, the variance coefficients (T) were established. 
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TABLE VIII.- CORRELATION COEFFICIENTS - R 
VARIANCE COEFFICIENTS - T 


ELEMENT 

SC 

oc 

T 

R 

R3 

X 


18.16555 

. 961076 

R5 

X 


18.13207 

.961537 

R1 

X 


17.90326 

.952979 

R6 

X 


17.03180 

.943269 

R7 

X 


16.96023 

. 940948 

L7 


X 

17 .98819 

.794069 

L6 


X 

17.87923 

.741361 

R7 


X 

17.98819 

.794069 

R6 


X 

17.87923 

.741361 

C2 


X 

17.32274 

.948035 

C7 


X 

15.67600 

.938431 

C6 


X 

15.44000 

.929121 

L7 

X 


17.28895 

.897051 

L6 

X 


17.11118 

.889213 

C7 

X 


14.62646 

.905976 

| C6 

X 


14.53590 

.897955 

R1 


X 

0 

0 

R3 


X 

0 

0 

C2 

X 


0 

0 

R5 


X 

0 

0 j 


Now to move away from the experimental situation, consider 
an actual network and assume that this network has undergone 
all the processes in items 1 through 5; that is, it has been 
simulated experimentally. 

Therefore, the desired response is known and any deviations 
from this response can be considered as failures. If a deviation 
is observed, it is desirable to establish the cause of these 
variances or in other words to isolate the failure. 
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The isolation procedure is as follows: 

1) F(0), the desired network response, is known. 







2 ) 

3) 


F (X) is the observed variant response. 


F (X) is compared to F(0) and, in the same manner as 
previously defined, a correlation coefficient, R* , and 
a variance coefficient, T* , are evaluated. 

4) Now consider the variance coefficients, T (k) , found 
experimentally from the modified transfer functions 
F(k). There exists a T for all possible single failures. 

5) Next, all variance coefficients, T(k), which are a 
certain specified percentage less than, and a certain 
specified percentage greater than, T* are designated. 

For example, the set (T(l), T(2), T(3), T(4), T(5), 

T ( 6 ) } when 


T 3 < T* £ T 4 . 

6 ) The set of variance coefficients, {T(l), T(2), T(3), 
T(4), T(5), T( 6 )} and also the variance coefficient T* 
have a corresponding correlation coefficient R* . 


7) Solve for the absolute value of the differences observed 
in the correlation coefficients (R(l), R(2), R(3), R(4), 
R(5), R ( 6 ) } of the modified transfer functions F(k): 
k = 1,2 ,... 6 and the correlation coefficient R* evaluated 
for the variant response F (X) . 


|R* - R(l) | = Z x 
| R* - R(2) | = Z 2 
|R* - R (3 ) | = Z 3 
| R* - R (4 ) | = Z 4 
|R* - R ( 5 ) I = Z c 

D 

|R* - R (6 ) I = Z c 
1 1 6 


8 ) The order of the numbers {Z]_, Z 2 , Z 3 , Z 4 , Z 5 , Zg} serves 
as the indicator of the likelihood of failure. The 
smallest Z has the highest probability of being the 
failure and so on. 


9) The entire process described above is applicable to the 
multi-failure situation and the algorithms for multi- 
failures will be considered later in this section. 
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Examples of the Isolation Process 

To illustrate the detection techniques we compile the 
expected response from the component data for the network of 
Figure 2 and examine the circuit after failure occurs. Two 
cases will be considered (Tables X and XI). 


TABLE IX.- FAILURE CASES 


W (i) 

BEFORE FAILURE 
F (0 , i) 

AFTER FAILURE 
CASE I 
F (X , i ) 

AFTER FAILURE 
CASE II 
F (X,i) 

10 1 - 5 

.362874 E-03 

.632455 E-08 

.156740 E-03 

io 2 

.362874 E-03 

.200000 E-07 

.156740 E-03 

io 2,5 

.362873 E-03 

.632455 E-07 

.156740 E-03 

10 3 

.362879 E-03 

.200000 E-06 

.156739 E-03 

IO 3 - 5 

.362824 E-03 

.632455 E-07 

.156736 E-03 

io 4 

.362372 E-03 

.199998 E-05 

.156705 E-03 

io 4 ’ 5 

.357991 E-03 

.632382 E-05 

.156393 E-03 

io 5 

.320148 E-03 

.199782 E-04 

.153314 E-03 

10 5 . 5 

.160346 E-03 

.629477 E-04 

.125939 E-03 

10 6 

.124238 E-04 

.342242 E-04 

.247579 E-04 


CASE I ; 

T* = 17.98726, R* = .794069 
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TABLE X.- DATA FOR ISOLATION PROCESS 


Case I 


SET A 

SET B 

SET C 


Corresponding 

Associated 

T* Interval 

Component 

R-Score 

16.96023 

R7-SC 

. 940948 

17.0313 

R6-SC 

. 943269 

17.11118 

L6-SC 

. 889213 

17.28895 

L7-SC 

. 897051 

17.32274 

C2-OC 

. 948035 

17.87923 

R6-OC 

.741361 

17 . 87923 

L6-OC 

.741361 

17.90236 

Rl-SC 

.952979 

17 . 98726 

R7-OC 

.794069 

17.98819 

L7-OC 

.794069 

18.13207 

R5-SC 

.961537 

18.16555 

R3-SC 

. 961076 


TABLE XI.- DATA FOR ISOLATION PROCESS - CASE II 


T = 17.32274 R = .948035 
2% Interval Specification for T gives 
LOWER LIMIT = 16.97629 
UPPER LIMIT = 17.66919 


SET A 

SET B 

SET C 


Corresponding 

Associated 

T* Interval 

Component 

R-Scores 

17.0317 

R6-SC 

.943269 

17.1118 

L6-SC 

.889213 

17.28895 

L7-SC 

. 897051 

17.32274 

C2-OC 

. 948035 
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To obtain a set of T scores all experimental variance 
coefficients within 6% (above or below) T* will be taken. It 
should be noted that T* will not always equal some T (k) score 
and R* will not always equal some R(k) score. 


The 6% specification gives a variance coefficient interval 
around T* with a lower limit of 16.06650 and an upper limit of 
19.0665. Now, from the experimental scores in Table VIII, the 
necessary data for fault isolation are illustrated in Table X. 
The following ordered interval around T* is obtained. 


Next, taking R* and Set C, the computation of | R* - R | = Z 
is performed 


| . 794069 
| .794069 
| .794069 
| .794069 
| .794069 
| .794069 
| . 794069 
| .794069 
j .794069 
| .794069 
| .794069 
I .794069 


. 940948 | 

= 

.146879 

.943269 j 

= 

.1492 

. 889213 | 

= 

.09514 

. 897051 J 

= 

.10298 

. 948035 | 

- 

.153966 

.74136l| 

= 

.052708 

.741361) 

= 

.052708 

.952979 | 

= 

.15891 

.794069 ) 

= 

0 

.794069 | 

= 

0 

. 961537 | 

= 

.167468 

. 961076 ) 

= 

.167007 


Z 

Z 

Z 

Z 

Z 

Z 

Z 

Z 

Z 

Z 

Z 

Z 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 


The order of parameter failure going from high probability 
to low is: { Z 9 = Zio, Zg = Z 7 , Z 3 , Z 4 , Zp, Z 3 , Z 5 , Zq, Zp2, Zip} 

which corresponds to {L7-OC = R7-OC , L6-OC = R6-OC, L6-SC , L7-SC, 
R7-SC, R6 -SC, C2-OC, Rl-SC, R3-SC, R5-SC} = the probability order 
of failure. The 6% interval for T* was used for illustrative 
purposes. Usually a smaller specification would be employed. 


I R* - R | = Z 

| .948035 - .940948) 
| .940835 - .889213 | 
| .948035 - .897051 | 
| .948035 - .948035) 


.007087 

.058822 

.050984 

.0000 
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Order of Parameter Failure by Probability 


= U 4 , z lr z 3 , z 2 } 

= { C2-0C , R 6 -SC , L7-SC , L 6 -SC } 


Simultaneous Failure of Several Elements 

For a single failure, the topology equation was expanded in 
two variables and H(P,Q) was solved for P in terms of element Q. 

For two simultaneous failures the expansion of H(P,Q,R) in 
three variables is employed where Q and R are network elements 
and P is the performance parameter. 

Solving H(P,Q,R) = 0 for 1/P defined to be F(0) yields 

1 _ H (P ,Q,R) + QH(P,Q',R) + RH(P,Q,R') + QRH(P,Q',R') 

P HCP'/^Rr + QH(P,Q',R) 4 - RH (P 1 , Q , R ' ) + QRH ( P ' , Q ' , R 1 ) 

When B parameters can be short-circuited or open-circuited 
simultaneously, 2B distinct modes of failure exist. Each mode 
depends on the combination of short-circuit and open-circuit. 

For B = 2 four formulas are listed. 


(1) 

H (P , Q , R) 

H (P"'7^,K) 

(3) 

H (P , Q 1 , R) 

H (P ' , Q ' ,K) 

(2) 

H (P , Q , R 1 ) 

(4) 

H (P , Q 1 , R ' ) 

H (P 1 ,U,R' ) 

H (P 1 , Q ' , R ' ) 


Since B = 2 elements, each Q parameter and each R parameter 
can function as either an impedance or as an admittance. 

There will be ( 2 ^ 6 ) combinations to be considered in 
Table XII. In general, the modified transfer functions can now 
be defined as F(k): k = 1 , 2 , 3 , . . . , 2 (^) where N is the number of 
elements that can fail and B is the number that can fail simul- 
taneously . 

The computations of the correlation coefficient and the 
variance coefficient do not change in the multiple failure 
situation. These coefficients are complied and utilized in the 
same manner as the single failure case. 
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TABLE XII.- POSSIBLE DOUBLE FAILURE CASES 



III.- FUTURE TASKS 

The technique presented may be extended as follows: 

a) Expand the computer program for multiple failures. 

Three simultaneous failures are a feasible goal. 

b) Test the technique for a wider range of circuits to 
establish its limitations and potential. 

c) Establish whether the coefficients can be utilized more 
effectively, particularly for multiple failures. 

The increasing number of components in microcircuits has 
established the necessity of failure analysis. The results of 
this analysis have provided one method for the isolation of 
failure in the frequency-dependent circuit. 
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